April 2003 
OHSTPY-HEP-T-03-005 



Phenomenology of the Minimal SO(10) SUSY Model 1 



Stuart Raby 

Department of Physics, The Ohio State University, 
174 W. 18th Ave., Columbus, OH 43210; 
On leave of absence, School of Natural Sciences, 
Institute for Advanced Study, Princeton, NJ 08450 



Abstract 

In this talk I define what I call the minimal SO(10) SUSY model. I then 
discuss the phenomenological consequences of this theory, vis a vis gauge and 
Yukawa coupling unification, Higgs and super-particle masses, the anomalous 
magnetic moment of the muon, the decay B s — > fi + fi~ and dark matter. 



1 Invited talk at PASCOS'03, Mumbai, India. January 3 - 8, 2003. This talk is based on work in 
collaboration with T. Blazck, R. Derrmsek, L. Roszkowski, R. Ruiz dc Austri and K. Tobe 



1 Minimal SO(10) SUSY Model 



Let me first define the minimal SO(10) SUSY model [MSO10SM] [1] and then I will 
discuss the phenomenological consequences of this theory. In the MSO10SM the quarks 
and leptons of one family are contained in a 16 dimensional spinor representation and 
the two Higgs doublets of the MSSM come from a single 10 dimensional representation. 
We have 




q= : , l= : , u 




16 D 



10 H D [H u , H d , T, T] 

For the third generation, there is a unique Yukawa coupling to the Higgs doublets with 

W D A 16 3 10 H 16 3 . 

As a consequence, the top, bottom, tau and v T Yukawa couplings satisfy Yukawa unifi- 
cation with X t = At = A T = Ap T = A. Note with a large Majorana mass for v T 
we have a see-saw mechanism resulting in a light left-handed neutrino, i.e. M p v T v T ==>- 
m Vr ~ mf/Mp. Although I will not discuss Yukawa terms for the first and second genera- 
tion of quarks and leptons, it is well known that it is not phenomenologically acceptable 
for them to receive all their mass via renormalizeable interactions with a single 10h- Nev- 
ertheless with effective higher dimensional interactions it is not difficult to obtain realistic 
fermion masses and mixing angles for all quarks and leptons [2] . Moreover if these mass 
matrices are hierarchical, we do not significantly affect the results derived from assuming 
exact Yukawa unification for the third generation. 

Finally, the soft SUSY breaking parameters are given by —C so f t = m\ & £f =1 16*16; + 
m 2 10 10* H 10 H -A A I63IO//I63 +M1/2 E? =1 (xiXi)+^B H u H d . All but one of these terms, 
are the most general consistent with SO(10). A universal scalar mass tuiq for all three 
families is an additional assumption. Hence, the soft SUSY breaking parameters are given 
by 

m 16 , m 10 , A , M1/2, tan p. 

Before continuing we note that one additional soft SUSY breaking parameter is needed, 
which we discuss next. 



1.1 Radiative EWSB with large tan/3 needs m 2 H < m 2 Hd 

It has been shown that there are two consequences of splitting the two Higgs doublet 
masses. It reduces the amount of fine tuning for radiative electroweak symmetry breaking 
[EWSB] [3]. In addition, it permits EWSB in an entirely new region of SUSY parameter 
space with m 16 ^> M1/2 [4]. 

We have considered the possibility of both D x term and "Just So" splitting [1]. In 
the former, we assume a soft SUSY breaking D term where Dx is the auxiliary field of a 
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U(l)x gauge interaction defined by SO (10) —> SU(5) x U(l)x- We then obtain 

m2 (H u ,H d )= rn 2 w T2D x 
m \q, u, e) = m w + D x 
m \d,L)= m 2 w -3D x . 

These boundary conditions at the GUT scale generically give the low energy result m~ < 
m~ which is bad for Yukawa unification. 
With "Just So" splitting we have 

™(H U , H d ) 
m (Q, u, e) 
m ll L) 

These boundary conditions give m| << m~ which is good for Yukawa unification. This 
latter case may be motivated by the fact that the Higgs multiplets must be special. They 
necessarily have a fi term and they also require doublet-triplet splitting. Moreover they 
have large threshold corrections at the GUT scale due to the tau neutrino. 

The v T contribution to Higgs splitting results from the Yukawa term (\ Vt v t L H u ) with 
\ Vt = X t = \ b = X T = A. Since v T couples only to H u , this contribution distinguishes 
H u and H d . At one loop we find Am 2 Hu w ^ (2m 2 w + m 2 + ^)log(^) + 

Taking typical GUT values for the parameters A = 0.7, M Vt = 10 14 GeV (which gives 
(Aml) atm ~ 10~ 2 eV 2 ), M G = 3 x 10 16 GeV and A 2 w 2m 2 w 4m 2 6 we obtain Am 2 H 
= ~ Arrijj u /m± ~ .07. This is "Just So" splitting of about the right size. 



m 
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1.2 Gauge coupling unification 



Presently, gauge coupling unification provides the only evidence for low energy SUSY [5, 
6, 7). 



a 2 




a--. 



Mi 



G 



Note, when threshold corrections are included, the three gauge couplings i = 1,2,3 
do not precisely meet at the GUT scale. Moreover, for consistency, one loop threshold 
corrections need to be included when using two loop RG running from Mq —> Mz- At 
one loop there are significant GUT threshold corrections from Higgs and GUT breaking 
sectors. We now define the GUT scale as the point where ai(M G ) = a 2 (M G ) = olq. A 
good fit to low energy data then requires e 3 = ^ a ^ M G)- a G) ^ _4%_ 



2 SO (10) Yukawa unification 

Let us now consider the constraint on the soft SUSY breaking parameters resulting from 
Yukawa unification [8, 1, 9]. Note, the GUT threshold corrections to Yukawa unification 
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from gauge and Higgs loops is typically insignificant. Weak scale threshold corrections, on 
the other hand, are proportional to tan/3 and cannot be ignored [10, 11]. The dominant 
contributions are given by 5m b /m b = Amf + Am^ + Am b ° 9 + ■ ■ ■ where the first comes 
from a gluino-sbottom loop, the second from the chargino-stop loop and the third from 
finite wave function renormalization graphs. Note, in general we have Araf ~ -AmJ > 
for \i > [ our conventions ] . The first two contributions are tan j3 enhanced and can be 
~ 50 %, while the typical size of the log contribution is ~ +6 %. The contribution to 
the top quark mass is not tan j3 enhanced and although the contribution to the tau mass 
is; nevertheless it is small due to the smaller values of the relevant gauge and Yukawa 
couplings at the electroweak scale. Finally, good fits to top, bottom and tau masses 
require 5mb/m b < — 2 



2.1 Data favors fi > 

We now argue that two pieces of low energy data favor positive values of /J. The first is 
the rate for the process b — > 57 and the second is the anomalous magnetic moment of the 
muon. In the first case, the chargino term typically gives the dominant SUSY contribution 
and for \i > it has opposite sign to the standard model and charged Higgs contributions, 
thus reducing the branching ratio. This is desirable since the SM contribution by itself 
is a little too large. As a result, trying to fit the data with \x < is problematic. In 
the second case, the contribution to the muon anomalous magnetic moment due to new 
physics (beyond the standard model) is measured to be a^ EW x 10 10 = 33.9 (11.2) [e + e~ 
- based] or 16.7 (10.7) [r - based] [12]. There are two results depending on whether one 
uses e + e~ or r data to determine the hadronic contribution to the amplitude. Note, in 
either case the sign of a^ EW is positive. Moreover in SUSY this sign is directly correlated 
with the sign of \i [13]. Again favoring positive /x. Hence we consider only positive \i in 
our analysis. 



2.2 x 2 Analysis 



We have performed a x 2 analysis of the MSO10SM with 11 input parameters defined at 
the GUT scale and 9 low energy observables in our \ 2 function [1]. 

The 11 input parameters at Mq are [A, aa, Mq, £3; "^10, A , tan/3(Mz), Dx 
[ D term splitting] ( or Am 2 H [Just So Higgs splitting]), mi 6 , /i, Mi/ 2 ], where the last 
three parameters are fixed while we vary 8 parameters using the CERN package Minuit 
to minimize x 2 . The 9 observables (experimental/theoretical uncertainty) [X^ xp (cr^)] 



defining x 2 



(X? 



-X 



theory y2 



P 



NEW , M z , M w , M t = 174.3 (5.1), 



are given by a, a s (M z ) 

m b (m b ) = 4.20 (0.20), M T \. 



0.118 (0.002), 



2.3 Bottom Line 

The bottom line result of our analysis is that Yukawa unification is possible only in a 
narrow region of SUSY parameter space. The result is also easy to understand. Since 
for ji > and bm b jm b < —2% we need |Am^| > AroJ. However Amf ~ ^ tan/3, 
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AmJ « tan/5 and An\ g « |& log(^-) ~ 6 %. In order to enhance the chargino 

contribution, we can make the numerator larger by making A t large and negative. This is 
accomplished by making Aq at Mq large and negative, i.e. A t « -<=>- Aq << 0. This 
also has the effect of making the denominator for the chargino contribution smaller since 

the stop mass matrix is of the form ^ m ^ ^2 * J ■ As a consequence we naturally 

obtain rrtf << m- b \ enhancing the chargino, in comparison to the gluino, contribution. Of 
course in order not to have a negative stop mass squared we need to make m 16 large. As 
a result of the \ 2 analysis we find that good fits require A ~ —2 mi 6 , mm ~ s/2 mi 6 , 
m i6 > 2 TeV 3> u, M1/2, an d Am^ ~ 10 %. In Fig. 1 we show the \ 2 contours for two 
different values of m\Q as a function of u and M\/2- It is clear that x 2 improves as we 
increase tuiq. Note also that the dominant pull for \ 2 is d ue to the bottom quark mass 
corrections as can be seen in Fig. 2. 
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Figure 1: x 2 contours for mi 6 = 1500 GeV (Left) and m\% = 2000 GeV (Right). The 
shaded region is excluded by the chargino mass limit m x + > 103 GeV. 



3 Summary - Minimal SO(10) SUSY Model 

Before discussing some phenomenological consequences of the MSO10SM, let us summarize 
the main ingredients of the model. We assume a supersymmetric SO (10) GUT with quarks 
and leptons in 16s. In addition we assume that the minimal Higgs content of the MSSM 
(H u , Hd) are contained in a single 10. Finally for the third family we assume the minimal 
Yukawa interaction with 

W D A 16 3 10 H 16 3 . 

The direct consequences of MSO10SM follow. 

• Gauge coupling unification - ag, Mq, €3 ~ —4%; 

• Yukawa unification - X t = A& = A T = Xp T = A; 
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Figure 2: Contours of constant ra&(rah)[GeV] (Left) and in % (Right) for m 16 = 2000 
GeV. 

• Soft SUSY breaking parameters 2 - m 16 , mio, ^4o; Myz, tan/3~50, Am 2 H ; 

• Satisfying - A ~ — 2 m 16 , m 10 ~ m 16 , mi 6 > 2 TeV ^> fi, M1/2, and 
Am^ ~ 10%. 

The last condition is required in order to fit the precision low energy electroweak data, 
including the top, bottom and tau masses. In addition to the above defining properties 
of the MSO10SM, we find two additional direct consequences of the model. The first is 
a "natural" inverted scalar mass hierarchy which ameliorates the SUSY flavor and CP 
problems. Secondly, the rates for proton decay due to dimension 5 operators are decreased. 
We discuss these two unexpected benefits below. 

3.1 Inverted Scalar Mass Hierarchy 

One way to ameliorate the SUSY flavor and CP problems is to demand that the 1st & 2nd 
generation squarks and sleptons are heavy with mass >> TeV, while the 3rd generation 
scalars are light with mass < TeV. This is easily seen by focusing on the most severe 
flavor and CP violating processes [14]. The best bounds are for processes involving the 
two lightest families. For example, we have [14] - 

• B(fi -> erj) < 1.2 x lO" 11 \($ 1 12 )ll\ < 2.1 x 10- 3 (m r (GeV)/100) 2 
or \(S[ 2 ) LL \ <0.8 (m r (TeV)/2) 2 ; 

• Am K < Exp. =>■ \/\Re(5f 2 ) 2 LL \ < 1-9 x 10^ 2 (m^(GeV)/500) 

or ^J\Re(5 d 12 ) 2 LL \ < 7.6 x lO" 2 K~(TeV)/2); 

2 In general mig would have a family index. We have made the additional assumption of a universal 
squark and slcpton mass. 
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• e K < Exp. \/\Im(5( 2 ) 2 LL \ < 1.5 x 10- 3 (m g -(GeV)/500) 

or y/\Irn(5f 2 ) 2 LL \ < 6.0 x 10~ 3 (m$(TeV)/2); 

• d e N ~ 2(100/m;-(GeV)) 2 sm$ AjjB x 10~ 23 e cm < 4.3 x 10~ 27 e cm sin^ AjB < 
4 x 1(T 4 x (m [ (GeV)/100) 2 

or sin<f> A ,B < 0.16 x (m f (TeV)/2) 2 . 

Although a significant degeneracy of the first and second generation squarks and sleptons 
is still required, it does not require serious fine tuning. In fact, the flavor and CP problems 
are now completely amenable to solutions using non-abelian family symmetries. The 
question one now faces is how to obtain an inverted scalar mass hierarchy with the ratio 
of scalar masses S satisfying S = fh\ 2 / m i ^ 1- 

One way of obtaining this inverted scalar mass hierarchy is to assume that it results 
from Planck/GUT scale physics. However, it was shown that an inverted scalar mass 
hierarchy can be generated "naturally" as a consequence of renormalization group run- 
ning [15]. This latter possibility requires specific soft SUSY breaking boundary conditions 
at Mq. In particular it was found that the following boundary conditions can lead to val- 
ues of S > 400 [15]. Surprisingly, these boundary conditions are the same required by 
Yukawa unification. 

• wJq = mfj = m 2 D = m\ = m% = m1 6 ; 

• A t = A b = A T = A ; 

• Mi = M 2 = M 3 = M 1/2 ; 

• m Hu = m Hd ee m 10 ; and 

• Al = 2 m\ Q = 4 mf 6 with m 16 » 1 TeV. 



3.2 Suppressing proton decay 

Nucleon decay rates are significantly constrained by data from Super-Kamiokande [16]. 
In particular the decay mode p — > K + + v T , due to dimension 5 operators, is typically 
the dominant decay mode. In the large tan (3 regime the dominant Feynman diagram is 
given below. 

This one loop integral results in a loop factor characteristically of order 



Loop factor 



I67H m 



16 



Note the loop factor is minimized in the limit fi, M1/2 tuiq. This limit is once again 
consistent with Yukawa unification. Moreover it is only consistent with radiative EWSB 
with split H u , Hd masses. 
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4 Phenomenology 

Let us now consider some predictions of the MSOi SM. 

4.1 Light Higgs mass 

First consider the light Higgs mass. In the MSSM the light Higgs mass has an upper 
bound of order 130 GeV. This upper limit is achieved for large tan/3. Moreover the large 
radiative corrections to the Higgs mass are dominated by heavy stop masses. In our case 
we have tan/3 ~ 50, however we have relatively light stop, sbottom, and stau masses. As 
a result we find [1] 

m h = 114 ± 5 ± 3 GeV. 

In Fig. 3 we show the light Higgs mass contours as a function of /x, M1/2 for two 
values of mi 6 . For a more detailed analysis of the light Higgs mass prediction, see [1]. 

4.2 Muon anomalous magnetic moment 

The anomalous magnetic moment of the muon scales as (fi M1/2 tan/3)/raf 6 . Since we 
have mi6 > 2 TeV, we find [1] 

a SUSY < q x 1Q -10 

4.3 x° LSP - Dark Matter 

When m 16 is large, the standard neutralino annihilation channels via squark/slepton ex- 
change diagrams are severely suppressed. This typically leads to an excess cosmological 
abundance of neutralinos with Q x h 2 ^> 0.3. However, in the large tan/3 regime the neu- 
tralino annihilation channel (x°X° A hadrons ) is significantly enhanced. In fact, 
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Figure 3: Contours of constant [GeV] (solid lines) with \ 2 contours from Fig 1 (dotted 
lines) for m 16 = 1500 GeV (Left) and m 16 = 2000 GeV (Right). 

this annihilation channel is so effective that, on resonance, Q x h 2 <C 0.01. Thus we find, 
on the sides of the broad resonance peak, cosmological abundances consistent with dark 
matter observations [17]. In Fig. 4 we present an analysis of dark matter abundances 
in the MSOioSM. The green band is the region with acceptable values of Q x h 2 . Note 
that we have also included contours of constant branching ratio B(B S — > fi + fi~). This is 
important since this process is extremely sensitive to the value of the CP odd Higgs mass 
m A . 



4.4 Large tan/? and Quark Flavor Violation 

It has been shown that in the large tan/3 regime there are significant one loop SUSY 
threshold corrections to CKM matrix elements [11]. Once these corrections are included 
in an effective two Higgs doublet model below the SUSY breaking scale, the Higgs cou- 
plings are no longer flavor diagonal [18, 19]. Hence the process B s — > can pro- 
ceed through s-channel CP odd Higgs exchange with a tan (3 enhanced branching ratio 
B(B S — > /i + /i~) oc tan/? 4 . The effective two Higgs doublet Yukawa coupling to down 
quarks is given below. The matrices A d *° 9 , AX l J (5Xd) are the zeroth order down quark 
Yukawa coupling in a diagonal basis and the one loop correction to the Higgs couplings 
due to gluino (chargino) loops. 



r ddH 



—djA A 



diag 



dm H. 



0* 



0* 



-d Li AA^ d Rj 
-d u 5\ l id R3 Hl + Yi.c. 

As a result of the chargino loop correction, which is proportional to the square of the up 
quark Yukawa matrix, we must re-diagonalize the down quark mass matrix 

v cos (3 



rn 



Diagonal 



V 



ydiag 



A\ d + SX d tanP 



V. 



fit 



V2 
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For large values of tan/3, this results in a significant correction to the CKM matrix [11]. 
We then obtain the following couplings to the neutral Higgs mass eigenstates h, H, A 
given by [18, 19]. 



*-FV 



U'. 

LJ/. 



Ft, I 
FSl 



R 



Pt 



F A P 



R 



+ F£* R 



j i 



d'H 
d'A, 



where 



~ 8X1(1 + tan 2 13) cos (3 cos(a - /3), 
~ 5A^'(1 + tan 2 /3)cos/3 sin(a-/3), 
F$ ~ 5X^(1 + tan 2 /3) cos j3. 

It is the flavor violating coupling which gives the direct B s A coupling [19]. Note, the 
branching ratio B(B S — > in the cosmologically allowed region is close to the CDF 

bound (see Fig. 4). In [17] we show that the process B s — > fi + fi~ may soon be observed. 



m 16 = 3TeV, m A =300 GeV 




M 1/s . (GeV) 



Figure 4: Contours of constant y 2 for m 16 = 3 TeV and = 300 GeV. The red regions 
are excluded by m x + < 104 GeV (below and to the left of a black solid curve), rrih < 111 
GeV (on the right) and by Q x h 2 > 0.3. To the right of the black broken line one has 
m>h < 114.4 GeV. The green band corresponds to the preferred range 0.1 < Q, x h < 0.2, 
while the white regions below (above) it correspond to Q x h 2 < 0.1 (0.2 < Q x h 2 < 0.3). 
Also marked are contours of constant BR(5 S — > fi + fi~). The blue regions in the lower 
two panels are excluded by BR(B S — > /z + /x~) > 2.6 x 10 -6 . 
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5 Two loose ends 



5.1 Fine tuning ? 

We have been considering large squark and slepton masses with tuiq > 2 TeV. Since we 
have a "natural" inverted scalar mass hierarchy, the third generation squarks and sleptons 
are typically lighter than a TeV. Hence the radiative corrections to the Higgs mass, in the 
effective low energy theory, are not large. For example the radiative corrections at the 
electroweak scale are of order 5m\ oc m \ and they are safe for m f < 1 TeV. 

However, there is still the question of whether radiative EWSB requires significant fine 
tuning. It has been shown that the Z mass is most sensitive to the value of the gluino 
mass, M 3 [20]. For example, with tan/3 = 35, the following relation was obtained [20] 

Mf = -1.5 11 2 + 5.0 M 2 /2 + 0.2 A\ + 1.5 m\ % - 1.2 m 2 Hu - 0.08 m 2 Hd + 

BUT, recall we have /x, Mi/ 2 <C ruie- Thus this problem is ameliorated somewhat, al- 
though it is not completely eliminated. (See also the talk by S. Pokorski, this conference.) 

Finally, as discussed earlier the fine tuning for EWSB in the regime of large tan (3 is 
of order l/tan/3 when one has Higgs mass splitting [3]. 

5.2 SUSY Breaking Mechanism ? 

We have found that Yukawa unification in the MSO10SM is only consistent with the low 
energy data in a narrow region of soft SUSY breaking parameter space. It is clear that 
this idea would be considerably strengthened if there was a mechanism which "naturally" 
broke SUSY in this way. Unfortunately, this is not the case for the known SUSY breaking 
mechanisms. 

For example, gauge mediated SUSY breaking [GMSB] gives A = at the messenger 
scale. This is bad. Yukawa deflected GMSB can have non-zero A$ proportional to Yukawa 
couplings. Perhaps this might work, however the standard gauge contribution would have 
to be suppressed. Moduli-dominated string SUSY breaking may be possible [15]. The 
generic formula for stringy SUSY breaking is given by [21] 

m 2 a = (1 + 3f£ • ml/2 

M1/2 ~ 
A af3l = ±y/3[l +11^ + ^ + 11^- Y^J 0* m 3/2 

where n a is the modular weight of the field a; Tj, (i — 1, ■ • • , 6) are moduli; 0j parametrize 
the direction of the goldstino in the Tj field space, and Y*^ = 2(ReTi) c^ln (h) a p where 
h a/ 3^ are the dimensionless Yukawa couplings. With the following modular weights [15]: 
n Q ,u,D,L,E,N = (0,-1/2,-1/2,0,0,0); n Hu>Hd = (-1/2,-1/2,0,0,0,0), and assuming 
(9 2 )* = (0, 0, 1/3, 2/3, 0, 0), and Y^ ~ 0, one finds m\ /2 = 2m 2 6 and 

A Q = ±2mi6, m io = V2m 16 . 
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The only problem with this idea is the absence of any existing string model with these 
properties. Thus the problem of a "natural" SUSY breaking mechanism consistent with 
MSO10SM Yukawa unification is the most urgent open theoretical question requiring fur- 
ther work. 

6 Summary 

In this talk I have defined the minimal SO(10) SUSY model and discussed some of its 
phenomenological consequences. The model predicts: 

• Gauge coupling unification with a G , M G , e 3 ~ —4%; 

• Yukawa unification with X t = A& = A T = \ Dt = A, and 

• Soft SUSY breaking parameters given by 3 m 16 , m 10 , A , M 1 / 2 , tan (3, Am 2 H . 

As a result of a \ 2 analysis [1] we find that the low energy precision electroweak data, 
including the top, bottom and tau masses, only gives good fits for soft SUSY breaking 
parameters satisfying: 

• A ~ — 2 m 16 , m 10 ~ y/2 m 16 , m 16 > 2 TeV ^> /i, Mi/ 2 , and Am 2 H ~ 10%. 
This region of parameter space has the virtue of giving: 

• a "natural" inverted scalar mass hierarchy which ameliorates the SUSY flavor and 
CP problems, and in addition 

• suppresses proton decay via dimension 5 operators. 

The MSOioSM makes the following predictions: 

• It gives tan (3 ~ 50 and a light stop. As a consequence we find [1] rrih = 114 ± 5 ± 3 
GeV; 

• The decay B s — > fi + /j," is enhanced and may be observable in the near future [19]; 

• The SUSY contribution to the muon anomalous magnetic moment is suppressed 
with al USY < 6 x 10" 10 [1]; and 

• Finally, it gives cosmologically acceptable abundances of neutralino dark matter [17] 
(see Roszkowski, this conference). 
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